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Most  attention  to  improving  vacuum  thermionic  energy  conversion  device  (TEC)  technology  has 
been  on  improving  electron  emission  with  little  attention  to  collector  optimization.  A  model  was 
developed  to  characterize  the  output  characteristics  of  a  TEC  where  the  collector  features  negative 
electron  affinity.  According  to  the  model,  there  are  certain  conditions  for  which  the  space  charge 
limitation  can  be  reduced  or  eliminated.  The  model  is  applied  to  devices  comprised  materials 
reported  in  the  literature,  and  predictions  of  output  power  and  efficiency  are  made,  targeting  the 
sub-1000  K  hot-side  regime.  By  slightly  lowering  the  collector  barrier  height,  an  output  power  of 
around  1  kW,  at  >20%  efficiency  for  a  reasonably  sized  device  (~0.1m2  emission  area)  can  be 
achieved.  ©2013  Author(s).  All  article  content,  except  where  otherwise  noted,  is  licensed  under  a 
Creative  Commons  Attribution  3.0  Unported  License.  [http://dx.doi.org/10.1063/L4826202] 


I.  INTRODUCTION 

Producing  more  energy  at  high  efficiency  is  a  constant 
challenge  facing  humanity.  Thermionic  energy  conversion 
has  a  role  to  play  due  to  its  many  advantages:  the  design  is 
simple,  there  are  no  moving  mechanical  parts,  and  there  is 
no  efficiency-reducing  direct  conduction  of  heat  across  the 
device.1'2 

Despite  the  advantages  of  thermionic  engines,  there 
have  been  technical  challenges  that  have  prevented  adoption 
on  a  large  scale.  The  challenges  have  been  dealing  with  the 
high  temperatures  required  to  achieve  a  reasonably  good  effi¬ 
ciency  and  having  to  overcome  the  effects  of  a  space  charged 
limited  mode. 

The  first  challenge  is  to  create  a  thermionic  energy  con¬ 
version  device  (TEC)  which  operates  at  a  reasonably  low 
temperature,  but  still  outputs  acceptable  power  at  high  effi¬ 
ciency.  In  order  to  get  the  same  electron  emission  from  a  ma¬ 
terial  at  lower  temperature,  one  must  either  find  a  way  to 
lower  the  thermionic  barrier  of  that  material  or  one  has  to 
find  a  different  material  with  a  lower  barrier. 

Naively,  one  might  think  that  lowering  the  emission  bar¬ 
rier  of  the  emitter  at  all  costs  is  the  way  to  improve  electron 
emission  and  therefore  TEC  output  power  density.  It  is  true 
that  electron  emission  dramatically  increases  by  lowering  the 
barrier,  but  at  the  device  level  output  power  will  not  neces¬ 
sarily  increase.  One  reason  is  that  more  emitter  output  cur¬ 
rent  results  in  a  greater  negative  space  charge  effect. 
Therefore,  emitter  barrier  lowering  quickly  hits  a  point  of 
diminishing  returns  due  to  the  self-limiting  nature  of  the  neg¬ 
ative  space  charge  effect.  The  second  general  issue  stems 
from  the  fact  that  output  power  depends  on  the  product  of 
output  current  and  the  operating  voltage.  Negative  space 
charge  issues  aside,  the  optimal  operating  voltage  of  the  de¬ 
vice  is  nominally  the  difference  in  barrier  heights  of  the  two 
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electrodes.  A  strategy  of  lowering  the  emitter  barrier  without 
consideration  of  the  collector  will  also  quickly  reach  a  limit 
of  effectiveness  as  the  value  of  the  emitter  barrier  approaches 
that  of  the  collector. 

The  second  challenge  is  the  so-called  negative  space 
charge  effect.  The  output  current  of  a  TEC  is  limited  by  the 
negative  space  charge  effect  because  electrons  traversing  the 
interelectrode  space  create  a  negative  charge  barrier  which 
blocks  lower  energy  electrons  from  reaching  the  collector. 
Decreasing  the  interelectrode  spacing  is  one  approach  to 
space  charge  mitigation. 

Generally  speaking,  the  emitter  barrier  should  be  greater 
than  the  collector  barrier.  The  emitter  barrier  should  be  low 
and  its  Richardson’s  constant  should  be  high  such  that 
enough  current  is  available  to  reach  the  desired  output 
power.  The  electrodes  should  be  close  enough  together  to 
mitigate  the  negative  space  charge  effect,  bearing  in  mind 
that  decreased  spacing  increases  the  engineering  difficulty  of 
device  fabrication. 

Recent  work  has  focused  on  improving  the  emission 
characteristics  of  the  emitter  electrode.  Nemanich,  Koeck, 
and  collaborators  have  reported  favorable  emission  parame¬ 
ters  from  nitrogen’  and  phosphorus  doped  diamond.4  Further 
improvements  in  electron  emission  occur  when  a  sample  is 
irradiated  with  light  via  a  phenomenon  called  photon- 
enhanced  thermionic  emission  (PETE).  Schwede  et  al.  first 
observed  PETE  from  GaN  (Ref.  5)  and  Sun  et  al.  have  also 
observed  it  from  diamond.6 

Regarding  space  charge  mitigation,  Hatsopoulos  and 
Kaye  demonstrated  at  considerable  difficulty  that  a  close¬ 
spaced  configuration  mitigates  space  charge.7  More  recently, 
Lee  and  collaborators  have  applied  MEMS  fabrication  tech¬ 
niques  to  decrease  the  interelectrode  distance.8  Smith  et  al. 
have  predicted  that  an  emitter  which  features  a  negative  elec¬ 
tron  affinity  (NEA),  such  as  hydrogen  terminated  diamond, 
could  mitigate  or  eliminate  the  negative  space  charge 
effect.9'10 
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Improvement  of  the  emission  properties  is  welcome 
advances  in  the  technology,  but  many  materials  already  exist 
with  quite  good  emission  at  low  temperatures:  scandate  elec¬ 
trodes,  BaO,  CsO,  and  other  oxide  materials"  all  have 
exceptional  electron  emission  properties. 

Little  investigation  is  being  conducted  on  improving  the 
collector  electrode  technology.  In  this  paper,  a  model  is 
developed  which  considers  the  effect  of  a  negative  electron 
affinity  collector  on  electron  transport  through  a  vacuum 
thermionic  engine.  It  is  shown  that  the  negative  electron  af¬ 
finity  reduces,  or  in  some  cases  eliminates  the  negative  space 
charge  effect.  The  model  is  applied  to  a  TEC  which  com¬ 
prised  materials  reported  in  the  literature  and  could  conceiv¬ 
ably  be  constructed  today.  Finally,  it  is  shown  that  of  all  the 
device  parameters,  the  output  power  density  and  efficiency 
can  be  most  effectively  optimized  by  lowering  the  collector 
barrier.  The  model  predicts  that  a  reasonably  sized  device 
(<lm2  emission  area)  can  produce  around  1  kW  at  >20% 
efficiency  at  a  hot-side  temperature  of  sslOOOK. 

II.  ELECTRON  TRANSPORT 

A  thermionic  engine  (aka  TEC)  is  a  vacuum  device  that 
converts  heat  directly  to  electrical  work.  It  is  a  heat  engine 
and  can  be  analyzed  as  such.  A  schematic  diagram  is  shown 
in  Fig.  1. 

Electron  transport  is  most  easily  understood  by  use  of 
the  electron  motive  and  electron  motive  diagram;  the  motive 
diagram  is  similar  to  the  band  diagram  in  a  solid-state  de¬ 
vice.  Examples  of  a  motive  diagram  in  the  case  where  the 
emitter  electrode  has  positive  electron  affinity  (PEA)  and  the 
collector  has  NEA  are  given  in  Fig.  2.  Subscripts  E  and  C 
denote  emitter  and  collector,  respectively.  The  electrodes  are 
separated  by  a  distance,  d.  The  motive  at  any  point  is 
denoted  by  i/r,  and  the  vacuum  level  at  either  surface  is 
denoted  by  ip  appropriately  subscripted  (e.g.,  in  the  case  of 
the  emitter:  i l/E).  The  maximum  motive  is  denoted  by  i//m. 
The  quantity  eV  represents  the  potential  of  the  collector  rela¬ 
tive  to  the  emitter,  V  is  referred  to  as  the  output  voltage  and 
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FIG.  1.  General  schematic  of  a  vacuum  TEC.  The  emitter  and  collector  elec¬ 
trodes  are  enclosed  in  a  vacuum  container  and  separated  by  some  distance. 
The  emitter  is  in  thermal  contact  with  a  thermal  reservoir  at  a  higher  temper¬ 
ature,  and  the  collector  is  in  thermal  contact  with  a  thermal  reservoir  at  a 
lower  temperature.  Electrons  are  thermionically  emitted  from  the  emitter, 
travel  across  the  interelectrode  space,  and  are  absorbed  by  the  collector.  The 
electrons  travel  from  the  collector  through  an  electrical  lead,  through  an 
external  load  where  work  is  done,  and  back  to  the  emitter  to  complete  the 
circuit. 


1.  Saturation  2.  Space  charge  3.  Virtual 
point  limited  mode  critical  point 


FIG.  2.  Various  conditions  of  electron  motive  in  a  TEC  where  the  collector 
features  NEA.  The  emitter  is  depicted  on  the  left  on  each  subplot,  and  the 
collector  on  the  right.  The  maximum  motive  is  indicated  by  a  “+”  symbol. 
The  output  voltage,  V,  is  shown  only  for  the  virtual  critical  point  to  avoid 
cluttering  the  other  two  motive  diagrams.  Coincidentally,  in  the  case 
depicted,  the  saturation  point  occurs  at  an  output  voltage  of  nearly  zero. 


e  is  the  fundamental  electronic  charge.  The  output  voltage  is 
physically  determined  by  the  output  current  and  the  external 
electrical  load  of  the  device.  For  the  case  of  PEA,  cf>  denotes 
the  thermionic  barrier  (sometimes  referred  to  as  work  func¬ 
tion),  and  in  the  case  of  NEA  the  barrier  is  denoted  by  £ 
which  is  the  difference  between  the  conduction  band  mini¬ 
mum  and  the  Fermi  energy.  The  value  of  NEA  is  denoted  by 
X-  The  quantity  i Pc  cbm 's  the  conduction  band  minimum  of 
the  collector,  and  /(  denotes  the  Fermi  energy. 

Electrons  inside  the  emitter  with  energy  greater  than  ipE 
are  thermionically  emitted  into  the  vacuum.  The  current  den¬ 
sity  of  the  emitted  electrons  is  given  by  Richardson’s 
equation 

i  =  ^exp(-|),  (1) 

where  ./  is  the  emission  current  [Acm  2|,  A  is  Richardson’s 
constant,  |Acm  2K  2|,  T  is  the  temperature  [K],  (p  is  the  bar¬ 
rier  height  or  work  function  [eV],  and  k  is  Boltzmann’s  con¬ 
stant  [8.62  eV  K-1]. 

Since  the  device  is  evacuated,  electrons  traveling  across 
the  interelectrode  space  represent  a  net  negative  charge.  For 
some  values  of  voltage,  this  negative  space  charge  can  cause 
the  motive  in  the  interelectrode  space  to  be  greater  than  the 
emitter  barrier,  presenting  an  additional  space  charge  barrier 
to  the  thermionic  electrons.  This  negative  space  charge  effect 
decimates  the  output  current  for  most  devices  and  is  a  signifi¬ 
cant  challenge  in  creating  a  viable  thermionic  device. 

III.  THEORY 

Electron  transport  in  the  face  of  the  negative  space 
charge  effect  is  modeled  using  a  Vlaslov-Poisson  system, 
following  the  success  of  this  approach  in  the  past.12'13 
Consider  a  TEC  with  a  PEA  emitter  and  a  NEA  collector. 
We  would  expect  this  device  should  experience  the  negative 
space  charge  effect  under  most  typical  operating  conditions, 
but  we  assume  the  collector  is  sufficiently  cool  so  that  back 
emission  is  negligible.  Electrons  arriving  at  the  collector 
with  energy  greater  than  t \/c  CBM  are  assumed  to  be  absorbed 
by  the  collector. 
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As  the  output  voltage  changes,  this  device  will  pass 
through  several  unique  modes  of  electron  transport.  The 
accelerating  mode  is  the  condition  where  the  maximum 
motive  occurs  at  the  emitter  electrode.  The  saturation  point 
is  the  greatest  value  of  output  voltage  such  that  the  maxi¬ 
mum  motive  occurs  just  outside  the  emitter;  it  represents  a 
bound  on  the  accelerating  mode  and  is  depicted  in  Fig.  2.1. 

As  the  output  voltage  increases,  the  maximum  motive 
exists  somewhere  within  the  interelectrode  space.  The  maxi¬ 
mum  motive  barrier  reduces  the  current  traveling  across  the 
device,  and  so  this  set  of  voltages  is  known  as  the  space 
charge  limited  mode.  The  motive  diagram  for  a  typical  point 
in  the  space  charge  limited  mode  is  depicted  in  Fig.  2.2.  The 
virtual  critical  point  occurs  when  the  maximum  motive  and 
collector  conduction  band  minimum  coincide  at  the  same 
height.  This  condition  is  shown  in  Fig.  2.3. 

As  the  voltage  increases,  the  conduction  band  minimum 
of  the  collector  limits  the  electron  current  entering  the  col¬ 
lector:  Electrons  with  energy  less  than  the  value  of  i/r c  cbm 
cannot  enter  the  collector  and  are  scattered  back  to  the 
emitter. 

At  some  value  of  voltage,  the  maximum  motive  occurs 
immediately  outside  the  collector.  This  situation  is  referred 
to  as  the  critical  point.  For  voltages  greater  than  the  critical 
point  voltage,  the  TEC  is  in  the  retarding  mode. 

Between  the  critical  point  and  the  virtual  saturation 
point,  the  analysis  of  Langmuir1-  exactly  models  the  electron 
transport.  For  output  voltage  above  the  virtual  saturation 
point,  the  electron  transport  differs  from  Langmuir  but  the 
current  is  ultimately  limited  by  the  position  of  t Ac  cbm- 
Therefore,  the  output  current  is  expressed  by  Eq.  (2)  which 
is  unaffected  by  the  precise  details  of  the  motive  in  the  inter¬ 
electrode  space 
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This  value  corresponds  to  a  value  of  ycs  according  to 
Langmuir’s  solution  to  the  dimensionless  Poisson’s  equation. 
From  Fig.  2. 1  and  the  definition  of  y 

eVs  =  (j^E  +  Xc  ~  Cc  ~  ycskTE-  (4) 


B.  Virtual  critical  point 

By  the  definition  of  the  virtual  critical  point,  and  refer¬ 
encing  Fig.  2.3 

Xc  =  VLj  —  •A CVR ■  (5) 

The  subscript  CVR  indicates  Collector  Virtual  critical  point. 
Using  the  definition  of  y  and  substituting 

Xcvr  =  T^r  ■  (6) 


Using  the  value  of  yCVR,  one  can  use  the  dimensionless  solu¬ 
tion  to  Poisson’s  equation  to  determine  £cvr-  From  the  defi¬ 
nition  of  ^ 


Zevr  —  £cvr  '  d 
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Using  the  dimensionless  Langmuir  solution,  yEVR  can  be 
determined  from  tEyR.  From  the  virtual  critical  point  motive 
diagram  and  the  definition  of  y 

eVVR  =  (j)E  -  Cc  +  kTEyEVR  ■  (8) 


/  =  ATE  exp  - 


Cc  +  eV\ 
kTE  ) 


(2) 


Notation  from  Langmuir’s  analysis  is  used  in  the  follow¬ 
ing  derivations  and  the  highlights  are  listed  here.  In  the  space 
charge  limited  mode,  Langmuir  converts  the  motive  and 
position  to  the  dimensionless  quantities;  y  =  is  the 

dimensionless  motive  and  £  =  (x  —  xm)  (2™tf  )  yy  is  the 

dimensionless  position.  The  quantity,  xm,  is  the  position 
where  i//,„  occurs.  It  is  worth  noting  that  the  appearance  of ./ 
in  the  expression  for  £  is  an  indication  that  we  are  self- 
consistently  solving  the  Vlaslov  and  Poisson  equations. 
These  dimensionless  quantities  allow  Langmuir  to  write 
Poisson’s  equation  in  a  universal  form  and  thus  numerically 
calculate  the  solution.  Hatsopoulos  and  Gyftopoulos  give  a 
clear  exposition  of  an  algorithm  to  determine  the  output  volt¬ 
age  given  a  value  of  space-charge  limited  current  density.14 
Thus,  the  output  current  density  vs.  output  voltage  can  be 
determined  for  a  set  of  operating  parameters. 


C.  Space  charge  limited  mode 

The  method  of  calculating  the  output  current  character¬ 
istic  in  the  space  charge  limited  mode  is  nearly  identical  to 
the  Langmuir  case,  one  difference  being  the  space  charge 
mode  extends  from  the  saturation  point  only  to  the  virtual 
critical  point  when  the  collector  exhibits  NEA.  The 
Hatsopoulous-Gyftopoulos  algorithm  to  find  the  output  volt¬ 
age  given  a  value  of  output  current  density  can  be  adapted  to 
this  case. 

1.  Given/,  calculate  yE  using  Eq.  (9). 

2.  Calculate  £E  using  Langmuir’s  solution  to  the  dimension¬ 
less  Poisson’s  equation. 

3.  Compute  £c  from  Eq.  (10). 

4.  Again  using  Langmuir’s  solution  to  the  dimensionless 
Poisson’s  equation,  calculate  yc  from  Qc. 

5.  The  output  voltage  is  given  by  Eq.  (1 1). 

7e  —  l'1  (7- ) ,  (9) 


A.  Saturation  point 

In  this  case  the  output  current  is  given  by  the  saturation 
current  of  the  emitter,  and  the  saturation  point  dimensionless 
distance  at  the  collector  is  given  by 


£c  —  £e  +  d 


/  2nme 


4k3 


1/4  ,1/2 

JVR 

I  7-3/4  ’ 


(10) 


eV  =  (pE  +  kTEyE  -  (£c  -  Xc  +  kTEyc). 


(ID 
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IV.  ANALYSIS 

A.  The  model  applied  to  a  plausible  TEC 

The  model  can  be  used  to  predict  the  performance  of  a 
TEC  created  from  materials  reported  in  the  literature. 
Consider  a  TEC  constructed  from  a  scandate  electrode 
described  by  Gartner  et  al. 1 1  as  the  emitter,  and  a  phospho¬ 
rous  doped  diamond  electrode  described  by  Koeck  et  al.4  as 
the  collector.  The  emitter  has  a  barrier  of  1.16eV  and  a 
Richardson’s  constant  of  7.8Acm  2K  2,  while  the  collector 
has  a  barrier  of  0.9  eV  and  a  Richardson’s  constant  of 
1  x  10  5Acm  2K  2.  Efficiency  of  the  device  can  be  calcu¬ 
lated  by  considering  the  electronic  heat  transport  and  the 
Stefan-Boltzmann  heat  transport;  unfortunately,  no  Stefan- 
Boltzmann  emissivity  data  are  quoted  for  either  material. 
The  scandate  material  is  usually  applied  as  a  thin  coating  to 
a  refractory  metal  such  as  tungsten,  and  the  diamond  is  CVD 
deposited  as  a  thin  film  on  molybdenum.  Thus,  a  mid-range 
emissivity  value  of  0.5  is  assumed  for  both  electrodes.  No 
value  is  quoted  for  the  NEA  of  the  phosphorus  doped  dia¬ 
mond  collector.  According  to  the  literature,  the  magnitude 
could  be  as  large  as  leV,15'16  so  a  conservative  value  of 
0.5  eV  is  chosen.  The  interelectrode  spacing  is  chosen  to  be 
10  /tin.  Finally,  assume  the  emitter  is  held  at  a  temperature 
of  1000K  while  the  collector  is  at  300  K.  At  such  low  tem¬ 
perature,  the  back  thermionic  emission  from  the  collector  is 
negligible,  consistent  with  the  model.  The  output  power 
characteristic  of  these  parameters  as  calculated  by  the  model 
is  given  by  the  dashed-dotted  line  in  Fig.  3.  The  maximum 
output  power  density  is  1.03Wcm~2  at  an  efficiency  of 
17.8%  and  occurring  at  an  output  voltage  of  0.38  V.  Thus, 
the  emission  area  of  the  device  would  be  0. 1  m2  to  produce 
1  kW  output  power.  The  output  power  density  and  emission 
area  of  this  device  are  reasonable,  but  the  efficiency  falls 
short  of  the  20%  target. 

In  the  physical  world,  materials  parameters  such  as 
Richardson’s  constant,  emissivity,  and  barrier  height  are  usu¬ 
ally  coupled.  The  computer  model  is  not  constrained  by  this 
coupling  and  so  individual  parameters  can  be  adjusted  to 
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FIG.  3.  Output  power  characteristic  showing  the  reduction  to  zero  of  the 
space  charge  limited  mode  by  decreasing  the  interelectrode  spacing.  Inset: 
Motive  diagram  at  the  condition  where  the  saturation  and  virtual  critical 
point  coincide. 
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FIG.  4.  Maximum  output  voltage  vs.  various  parameters  with  all  others  held 
fixed.  From  left  to  right:  emitter  temperature,  emitter  barrier  height,  collector 
barrier  height,  interelectrode  distance,  and  collector  NEA.  The  red  circles 
indicate  the  common  starting  set  of  parameters. 


quickly  see  which  parameters  have  the  biggest  effect  on  de¬ 
vice  performance.  This  simulation  is  valuable  to  build  intu¬ 
ition  about  how  this  system  behaves  since  the  system  does 
not  have  a  closed-fonn  solution.  Using  the  device  just 
described  as  a  starting  point,  individual  input  parameters  are 
adjusted  (keeping  all  others  fixed)  to  calculate  the  maximum 
efficiency.  Fig.  4  depicts  the  results  of  such  a  calculation  as 
an  array  of  plots.  The  plots  depict  the  maximum  efficiency 
vs.  a  particular  input  parameter;  the  plots  from  left  to  right 
show  the  effect  of  adjusting  the  emitter  temperature,  emitter 
barrier,  collector  barrier,  collector  NEA,  and  interelectrode 
spacing,  respectively.  The  red  circles  indicate  the  value  of 
the  common  starting  set  of  parameters. 

The  first  plot  shows  that  the  efficiency  increases  as  emit¬ 
ter  temperature  increases,  crossing  20%  efficiency  at 
Te=  1070  K.  This  result  is  not  surprising  because  increasing 
the  emitter  temperature  increases  the  number  of  electrons  in 
excited  states  and  thus  more  electrons  are  available  to  cross 
the  device. 

The  second  plot  shows  that  efficiency  increases  rapidly  as 
the  emitter  barrier  decreases,  then  levels  off.  If  the  emitter  bar¬ 
rier  is  high,  few  electrons  escape.  The  output  current  is  small, 
and  therefore  the  negative  space  charge  effect  is  negligible. 
Additionally,  since  the  output  current  is  small,  both  the  output 
power  and  efficiency  will  be  small  as  well.  As  the  emitter  bar¬ 
rier  decreases,  the  output  current  increases  but  so  does  the 
negative  space  charge  barrier.  At  some  point,  the  emitter  bar¬ 
rier  is  low  enough  that  the  current  is  limited  mostly  by  the 
space  charge  and  not  the  emitter  barrier;  further  lowering  the 
emitter  barrier  has  a  diminishing  effect  on  the  output  current. 
The  space  charge  limitation  on  current  limits  the  output  power 
and  thus  efficiency.  Lowering  the  emitter  barrier  is  akin  to 
integrating  over  a  larger  interval  of  the  thennal  distribution  of 
the  electrons  of  the  emitter,  in  contrast  to  increasing  the  emit¬ 
ter  temperature  which  increases  the  scale  of  the  distribution. 

The  third  plot  shows  that  the  maximum  efficiency 
increases  as  the  collector  barrier  is  lowered,  crossing  20%  at 
Cc  =  0.85eV.  This  condition  yields  an  output  power  density 
of  1.17Wcm  2  (corresponding  to  an  emission  area  of 
0.086  m2  for  1  kW  total  output)  at  an  output  voltage  of 
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0.43  V.  Examining  the  motive  diagram  is  the  easiest  way  to 
understand  the  relationship  between  collector  barrier  and 
maximum  efficiency.  Consider  that  the  efficiency  is  calcu¬ 
lated  as  the  output  power  divided  by  the  rate  of  heat  input, 
and  the  output  power  is  the  product  of  the  output  current  and 
voltage.  In  the  motive  diagrams  depicted  in  Fig.  2,  the 
motive  due  to  the  negative  space  charge  effect  is  calculated 
using  the  vacuum  levels  of  both  electrodes  as  boundary  con¬ 
ditions.  On  the  collector  side,  the  vacuum  level  is  determined 
relative  to  the  emitter  Fermi  energy  by  the  output  voltage 
(times  the  fundamental  charge),  the  collector  barrier,  and  the 
collector  NEA.  For  a  smaller  value  of  collector  barrier,  the 
system  would  have  the  same  boundary  conditions  if  the  out¬ 
put  voltage  was  increased  (assuming  constant  value  of 
NEA).  This  tradeoff  does  not  affect  the  output  current  of  the 
device,  but  the  output  power  would  increase  as  a  result  of  the 
increased  voltage.  Moreover,  the  efficiency  would  increase 
as  well.  This  relationship  between  collector  barrier  and  effi¬ 
ciency  is  not  unique  to  the  case  where  the  collector  features 
NEA — a  similar  relationship  holds  in  a  traditional  vacuum 
TEC.  Indeed,  the  performance  of  a  TEC  has  been  analyzed 
in  terms  of  the  so-called  barrier  index,17'18  Vb  =  Vd  +  4>c- 
The  quantity  Vrl  is  generally  understood  to  be  the  plasma  arc 
drop  in  a  vapor  device,  but  could  be  interpreted  as  the  space 
charge  barrier  in  a  vacuum  device.  This  analysis  essentially 
holds  in  the  case  where  the  collector  features  NEA;  i.e.,  the 
efficiency  can  still  be  improved  by  decreasing  the  collector 
barrier  in  the  expression  for  the  barrier  index.  The  difference 
is  that  the  presence  of  collector  NEA  adds  a  parameter  to  in¬ 
dependently  adjust  the  space  charge  barrier  component  of 
the  barrier  index.  According  to  the  theory  developed  in  this 
report,  a  TEC  with  collector  NEA  will  out-perform  a  TEC 
without  NEA  as  a  result  of  space  charge  mitigation. 

Strictly  speaking,  collector  barrier  reduction  has  a  limit  of 
effectiveness  because  at  some  point  the  collector  will  experi¬ 
ence  appreciable  back  emission.  Back  emission  is  negligible 
over  the  entire  range  of  collector  barrier  height  values  shown 
in  Fig.  4.  The  back  current  attains  a  value  of  1  x  10-7  A  at 
(c  =  0.65eV,  which  is  many  orders  of  magnitude  smaller  than 
the  output  current. 

The  fourth  and  fifth  plots  show  that  efficiency  tapers  off 
as  the  value  of  NEA  increases  and  the  interelectrode  distance 
decreases,  respectively.  These  results  are  two  manifestations 
of  the  same  phenomenon  having  to  do  with  space  charge  mit¬ 
igation  as  explained  in  Sec.  IV  B.  Leaving  the  details  for 
Sec.  IV  B,  when  the  collector  exhibits  NEA,  there  exists  a 
set  of  device  parameters  such  that  the  negative  space  charge 
effect  is  completely  eliminated.  Once  the  space  charge  is 
eliminated,  the  device  operates  at  its  ideal  performance; 
decreasing  the  interelectrode  distance  or  increasing  the  NEA 
has  no  further  effect  on  the  performance  of  the  device,  and 
therefore  the  efficiency  does  not  improve  further. 

B.  Reduction/elimination  of  the  space  charge  mode 
via  collector  NEA 

According  to  the  model,  there  are  some  cases  for  which 
a  TEC  employing  a  collector  featuring  NEA  experiences  no 
space  charge  limited  mode.  It  is  well  known  that  the  negative 


space  charge  effect  can  be  mitigated  by  reducing  the  intere¬ 
lectrode  spacing  of  a  vacuum  thermionic  engine;  space 
charge  mitigation  in  a  device  with  a  NEA  collector  is  similar 
to  space  charge  mitigation  in  a  non-NEA  device.  The  differ¬ 
ence  is  that  in  the  non-NEA  device,  full  space  charge  elimi¬ 
nation  occurs  in  the  limit  of  d  — >  0,  while  in  a  NEA  collector 
device  full  space  charge  elimination  can  occur  at  a  finite  dis¬ 
tance.  The  initial  device  considered  in  this  report  will  experi¬ 
ence  no  space  charge  limited  mode  at  an  interelectrode 
distance  of  3.89  qm  and  below.  Lee  et  al.  have  shown  for 
very  small  values  of  interelectrode  spacing  the  efficiency  of 
a  TEC  degrades  due  to  near- field  radiative  coupling.19 
According  to  the  model  described  in  this  paper,  a  TEC  fea¬ 
turing  a  negative  electron  affinity  collector  can  achieve  ideal 
operation  with  an  interelectrode  spacing  sufficiently  large  to 
avoid  near-field  effects. 

Before  delving  into  a  derivation  of  the  conditions  under 
which  a  TEC  with  a  NEA  collector  avoids  the  space  charge 
limited  mode,  a  qualitative  sense  of  the  phenomenon  can  be 
obtained  by  individually  considering  the  effect  of  both  yc 
and  d  have  on  the  motive. 

First,  consider  a  TEC  with  a  NEA  collector  operating  in 
the  space  charge  limited  mode;  this  TEC  will  have  a  motive 
diagram  similar  to  the  motive  diagram  shown  in  Fig.  2.2. 
Increasing  yc  is  equivalent  to  lowering  the  position  of  the 
collector  vacuum  level,  i \)c.  Since  t//c  sets  the  boundary  con¬ 
dition  of  Poisson’s  equation,  the  solution  of  Poisson’s  equa¬ 
tion,  i.e.,  the  motive,  must  also  decrease.  For  simplicity  of 
the  argument,  assume  the  collector  Fermi  level,  qc,  is  suffi¬ 
ciently  low:  At  some  point  /c  will  be  sufficiently  large  that 
the  maximum  motive  will  occur  just  outside  the  emitter.  At 
this  point,  there  is  no  longer  a  space  charge  barrier,  and  the 
negative  space  charge  effect  has  been  mitigated.  Further 
increasing  y(:  has  no  additional  effect  on  the  output  current 
density  since  the  output  current  density  is  already  at  satura¬ 
tion.  A  corollary  of  this  result  pertaining  to  the  device  effi¬ 
ciency  is  depicted  in  the  fifth  plot  of  Fig.  4. 

Next,  consider  what  happens  to  the  space  charge  lim¬ 
ited  motive  as  the  interelectrode  distance,  d,  is  adjusted. 
This  explanation  is  based  on  the  observation  that  the  gen¬ 
eral  shape  of  the  motive  curve  is  fixed  (for  output  voltages 
less  than  or  equal  to  the  virtual  critical  point  voltage).  The 
general  shape  of  the  motive  curve  is  fixed  because  the  solu¬ 
tion  to  Poisson’s  equation  can  be  written  in  a  dimensionless 
form  as  given  by  Langmuir.  This  dimensionless  solution  is 
transformed  to  real  space  via  scaling  by  parameters  specific 
to  the  case  under  consideration,  but  the  scaling  from  dimen¬ 
sionless  motive  and  position  to  real  motive  and  position  is 
the  same  regardless  of  the  voltage  or  interelectrode  distance 
as  can  be  seen  in  the  definition  of  y  and  £  given  in  Sec.  III. 
Indeed,  the  algorithm  for  determining  the  output  current 
density  vs.  output  voltage  given  in  Sec.  IIIC  as  well  as 
Hatsopoulos-Gyftopoulos  algorithm  is  based  on  the  fact 
that  the  general  shape  of  the  motive  curve  is  fixed.  These 
algorithms  start  by  determining  the  extremum  of  the  fixed 
motive  curve  in  real  space,  then  mapping  that  point  to 
dimensionless  space,  and  finally  finding  the  points  on  the 
fixed  motive  curve  that  correspond  to  the  boundary  condi¬ 
tions  in  real  space. 
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Equipped  with  the  knowledge  that  the  general  shape  of 
the  motive  curve  is  fixed,  consider  the  motive  diagram 
shown  in  Fig.  2.1.  If  the  interelectrode  distance  is  decreased, 
the  output  voltage  must  increase  in  order  for  the  boundary 
condition  of  the  collector  to  be  satisfied  since  the  shape  of 
the  motive  curve  is  fixed.  As  d  is  lowered  in  this  way,  for 
some  value  of  d  the  output  voltage  will  be  such  that  i//c/;M  is 
coincident  with  i j/E,  and  the  TEC  will  experience  no  space 
charge  limited  mode.  Put  another  way,  since  the  shape  of  the 
motive  curve  is  fixed,  as  d  decreases,  the  collector  attributes 
(i 1>cbmi  •Ac i  A <c>  etc-)  travel  up  along  the  motive  curve  such 
that  i jjc  always  meets  the  curve.  If  one  considers  the  case  of 
the  virtual  critical  point,  a  similar  argument  can  be  made. 

A  consequence  of  the  above  analysis  is  that,  as  the  intere¬ 
lectrode  distance  decreases,  the  values  of  the  saturation  point 
voltage  and  virtual  critical  point  voltage  approach  the  same 
value.  This  result  is  depicted  in  Fig.  3  where  the  output  power 
characteristic  is  calculated  at  different  values  of  interelectrode 
distance.  The  effect  of  yc  and  d  on  space  charge  mitigation  is 
different  manifestations  of  the  same  phenomenon  because  the 
general  shape  of  the  motive  curve  is  fixed.  Since  the  shape  is 
fixed,  the  algorithm  to  determine  output  current  density  vs. 
output  voltage  amounts  to  locating  the  point  of  maximum 
motive  and  the  boundary  conditions  on  the  fixed  curve  accord¬ 
ing  to  some  constraints.  Both  yc  and  d  are  built  into  those 
constraints,  but  not  into  the  transform  from  dimensionless 
space  to  real  space.  In  general,  decreasing  yc  reduces  the  size 
of  the  voltage  interval  over  which  the  TEC  experiences  the 
space  charge  limited  mode,  as  does  reducing  the  interelectrode 
distance.  Therefore,  a  TEC  could  be  designed  to  have  a  larger 
interelectrode  distance  and  not  suffer  a  space  charge  limited 
mode  so  long  as  the  value  of  y(:  was  sufficiently  large. 

C.  Derivation  of  space  charge  elimination  condition 

To  quantitatively  formalize  the  previous  exposition, 
there  exists  a  volume  in  parameter  space  such  that  for  points 
inside  the  volume  and  on  the  surface  bounding  the  volume, 
there  is  no  range  of  voltages  for  which  the  TEC  passes 
through  a  space  charge  limited  mode.  These  entities  are 
referred  to  as  the  ideal  volume  and  ideal  surface.  For  some 
value  of  d,  called  didea/,  the  critical  and  virtual  saturation 
points  of  the  TEC  coincide;  this  combination  of  parameters 
lies  on  the  aforementioned  ideal  surface  in  parameter  space. 
At  this  point,  the  negative  space  charge  effect  is  completely 
eliminated  and  there  is  no  range  of  voltages  for  which  the  de¬ 
vice  passes  through  a  space  charge  limited  mode. 

The  inset  of  Fig.  3  depicts  the  general  motive  diagram 
corresponding  to  a  point  on  the  ideal  surface.  To  derive  this 
condition,  consider  a  TEC  with  parameters  on  the  ideal  sur¬ 
face:  It  is  equivalent  to  say  the  saturation  point  and  virtual 
critical  point  of  that  TEC  coincide.  We  first  compute  the  quan¬ 
tities  at  the  saturation  point.  From  the  figure,  ycs  becomes 

<12) 

So  Eq.  (4)  reduces  to 

eVs  =  4>e-  Cc-  (13) 


Note  also  that  the  output  current  density  is  equal  to  the  emit¬ 
ter  saturation  current  density. 

Next,  we  compute  quantities  of  the  virtual  critical  point. 
From  the  figure,  it  is  clear  that  the  output  current  density 
must  also  be  equal  to  the  emitter  saturation  current  density. 
Moreover,  ^EVR  equals  zero  and  therefore  yEVR  does  as  well. 
Equation  (8)  becomes 


eVVR  =  <f>E  -  Cc-  (14) 


The  output  voltage  and  output  current  density  of  both  the  sat¬ 
uration  and  virtual  critical  points  are  identical  and  therefore 
coincide. 

Another  consequence  of  the  fact  that  CEVR  =  0  and 
JVr  =  Jes  is  that  Eq.  (7)  can  be  rewritten  as 


3/4  ■ 


(15) 


Note  that  the  left  hand  side  of  Eq.  (15)  depends  on  the 
value  of  yCVR  via  Langmuir’s  solution  to  the  dimensionless 
Poisson’s  equation.  If  the  above  condition  is  met,  the  satura¬ 
tion  point  coincides  with  the  critical  point  and  the  TEC  is  on 
the  ideal  surface  in  parameter  space. 

Notice  also  that  there  are  device  parameters  on  which 
the  left  hand  side  of  Eq.  (15)  depends,  but  for  which  the  right 
hand  side  does  not.  For  example,  yc  affects  the  value  of  £CVr 
on  the  left  hand  side  but  no  quantity  on  the  right  hand  side.  If 
a  parameter  is  changed  such  that 


icvR  >  d 


^2mnre2\  ^  JXJ2 


elk3 


ES 
3/4  ’ 


(16) 


then  the  device  is  operating  within  the  ideal  volume  where  it 
experiences  no  space  charge  limited  mode.  If  the  inequality 
is  reversed,  the  device  will  pass  through  a  space  charge  lim¬ 
ited  mode  over  some  range  of  voltages. 


V.  CONCLUSIONS 

A  model  was  developed  which  considers  the  space 
charge  limited  electron  transport  through  a  TEC  where  the 
collector  features  NEA.  Using  this  model,  it  is  shown  that 
NEA  can  mitigate,  and  in  some  cases  eliminate,  the  nega¬ 
tive  space  charge  effect.  Calculations  were  performed  using 
the  model  to  show  output  characteristics  of  devices  made 
from  available  materials.  Using  material  parameters  quoted 
in  the  literature,  and  making  reasonable  inferences  on  miss¬ 
ing  parameter  values,  the  output  power  and  efficiency  of  a 
plausible  device  were  calculated  to  be  1.03Wcm~2  and 
17.8%,  respectively.  Since  the  efficiency  of  the  device  was 
below  the  target  of  20%,  the  input  parameters  were  individ¬ 
ually  adjusted.  The  device  can  reach  the  20%  target  when 
the  collector  barrier  falls  below  0.85  eV.  At  that  point,  the 
output  power  density  is  1.17Wcm~2.  Such  a  device  would 
output  1  kW  with  0.086  m2  emission  area.  Increasing  the 
emitter  temperature  to  1070K  would  also  achieve  20% 
efficiency. 
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VI.  COLOPHON 

The  model  was  implemented  as  a  module  written  in 
python.  Revision  0.3.2,  commit  3441d96a  (Ref.  20)  of  the 
software  implementing  the  model  was  used  for  the  calcula¬ 
tions  in  this  manuscript  along  with  scipy  version  0.11.0, 
numpy  version  1.7.0,  matplotlib  version  1.2.0,  and  python 
version  2.7.3 — all  installed  using  the  Homebrew  OS  X  pack¬ 
age  manager.2"  The  code  was  executed  on  a  MacBook  Air 
model  A1466  running  OS  X  10.8.4  (12E55),  Darwin  kernel 
version  12.4.0  with  a  2  GHz  Intel  core  i7  and  8GB 
1600  MHz  DDR3  memory.  The  software  to  perform  the  cal¬ 
culations  in  this  paper  was  developed  according  to  the  best 
practices  advocated  by  the  Software  Carpentry  project.21 
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